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Abstract: We consider higher derivative gravity lagrangians in 3 and 4 dimensions, which 
admit simple c-theorems, including upto six derivative curvature invariants. Following a 
suggestion by Myers, these lagrangians are restricted such that the fluctuations around (anti) 
de Sitter spaces have second order linearized equations of motion. We study c-theorems both 
in the context of AdS/CFT and cosmology. In the context of cosmology, the monotonic 
function is the entropy defined on the apparent horizon through Wald's formula. Exact black 
hole solutions which are asymptotically (anti) de Sitter are presented. An interesting lower 
bound for entropy is found in de Sitter space. Some aspects of cosmology in both D = 3 and 
D = 4 are discussed. 
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1. Introduction 

It is usually difficult to deal with equations of motion with more than two derivatives. Typi- 
cally such theories are plagued with ghosts ||| . String theory has in principle, a systematic 
way to compute higher derivative corrections to the low energy effective action. For instance 
it is well known that heterotic string has a correction that is quadratic in curvature [||, H|. If 
one computed fluctuations with this effective action, one would find ghosts. A conservative 
viewpoint to this problem is to argue that field redefinitions allow us to write the quadratic 
correction as a Gauss-Bonnet term in which case this problem disappears j|] . It is difficult to 
imagine that something similar can be done order by order around an arbitrary background. 
Of course one could take the attitude that there are an infinite set of such terms and once 
all of them are taken into account the mass of the problematic ghost modes would be pushed 
to infinity or essentially they would be removed from the spectrum. This would only work if 
the higher derivative terms are treated perturbatively. 

In this paper, we will consider higher derivative lagrangians in D = 3 and D = 4 that 
include upto six derivative terms, i.e., the action schematically reads 

]D „ / -/ o 1 A i \ E>2 , ~ r>3 



Iocj d u x^f=g{R-2h + \W + {lR 6 ), (1.1) 

where R 2 and R 3 denote a general set of four and six derivative curvature corrections. The 
effect of stringy R 2 corrections were first studied in While R 2 corrections arise in super- 
symmetric string theory, i? 3 corrections arise in nonsupersymmetric string theories ||. Six 
derivative theories have featured recently in (?], ||, |9|, [l(J 11 ] . Following a suggestion by Myers 



[0, [3^] , one of the main goals of this paper is to construct general lagrangians such that when 
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one considers fluctuations around (anti) de Sitter space, the linearized equations of motion are 
two derivative. We will not treat the higher derivative terms perturbatively in A or jl. One of 
the main reasons for not treating the higher derivative corrections perturbatively is to leave 
the possibility of probing issues such as lowering the viscosity bound [l^, [R]] in consistent toy 
models open. We will also demand that the resulting lagrangians admit for simple c-theorems 
in AdS/CFT |Q| and cosmology. If there is a flow between two theories, then the value of 
the c-function at the fixed points is supposed to be a measure for the number of degrees of 
freedom at the fixed points. 

The degrees of freedom on the CFT side are captured by central charges. In 1+1 di- 
mensional CFTs, unitarity, a conserved stress energy tensor and the Euclidean group of 
symmetries are enough to show the existence of a c-theorem whereby the ultraviolet value 
of the central charge is greater than the infrared value Using the AdS/CFT correspon- 
dence, it is straightforward to show this c-theorem using gravity equations of motion and 



the null energy condition [14]. In [||, postulating the existence of a simple c-theorem was 



used to derive the new massive gravity model |Tf| [[(], 17, [L8|, 19, |2^, pH}] . These theorems 



were considered further in [22]. In 23], c-theorems in arbitrary dimensions were investigated 



and an entanglement entropy interpretation for the quantity that was flowing was given in 
arbitrary dimensions. Related issues have been further discussed in ][24||. In the derivation of 
the theorem, the equations of motion and null energy conditions are used-the cosmological 
constant does not enter in an important way. It is natural to ask if a similar theorem exists 
in the context of cosmology. 



A c-theorem for cosmology similar in spirit [ 14 1 to the AdS / CFT correspondence was first 



proposed by Strominger in [ 25 1 . Here the starting point is to assume that the bulk metric is 
given by 

ds 2 = -dt 2 + a(t) 2 dyL 2 D _ x , (1.2) 

such that the scale factor at early and late times behaves like ^ — > Hi, ^ — > Hf respectively, 
with Hi being the inflationary era Hubble constant and Hf being the current value of the 
Hubble constant. In the intermediate stages, the evolution is governed by standard FRW 
equations. In two derivative Einstein gravity, the FRW equations give us 

H= ?L-£)* = -8nG(p + P), (1.3) 
a a 

so that if the null energy condition p + P > holds then 

1 

will be an increasing function in time 1 . Here s(t) is interpreted as the entropy on the apparent 
horizon. We will be interested mostly in the flat FRW case where the apparent horizon is the 



«(*) ~ (1-4) 



1 In the context of de Sitter/CFT correspondence, late times in the bulk correspond to ultraviolet in the 
supposed boundary field theory while early times correspond to the infrared. Bulk time evolution is thought 
to be an inverse RG flow (IR to UV) and hence corresponds to "integrating in" degrees of freedom. Hence it 
is natural to expect that entropy increases with time. 
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Hubble horizon. When the space is de Sitter, the notion of cosmological horizon and apparent 
horizon coincide. Thermodynamics aspects of apparent horizons have been studied in 26] 
and extended to Gauss-Bonnet and Lovelock theories in higher dimensions. In the context of 
dS/CFT, c-theorems have been investigated in j^7|. Of course as it stands eq.(|Q|) does not 
prevent a(t) from running off to zero in the early past. It is natural to expect that quantum 
corrections will be very important in this case and there will be corrections to eq.( |1.4| ). 

A first step towards understanding quantum effects is to study the inclusion of higher 
derivative curvature corrections in the action. Typically this may lead to problems for instance 
with ghosts although it is possible that the approach of [U in dealing with ghosts will still 
allow us to extract useful physical information. In 2+1 dimensions, an interesting higher 
derivative gravity theory was proposed in [|T6| , 28, ^] where four derivative terms R a bR ab — 
3/8R 2 are added to the Einstein-Hilbert lagrangian. The propagating degree of freedom 
is a massive spin-2 field. Around flat space, this degree of freedom can be shown to be 
unitary in spite of the appearance that the equations of motion are intrinsically higher order. 
When the equations of motion are higher order, schematically the graviton propagator looks 
like l/(p 2 (p 2 + m 2 )) cx 1/p 2 — l/(p 2 + m?) so that there is always an additional degree of 
freedom which the wrong sign kinetic term. In 2+1 dimensions it is possible to make the 
non-propagating mode have the wrong sign while the propagating massive mode have the 
right sign [28]. While this seems to work around flat space, this construction is problematic 
in the context of AdS/CFT. In this case, demanding that the massive mode is unitary in the 
bulk leads to the dual CFT have negative central charge [19]! 

One possible way out of this is to add more terms to the lagrangian as in || and tune 
the coefficients such that the equations of motion for the fluctuations around AdS space is 
two derivative pC[| . Then for a specific choice of the parameters it is possible to show that 
the CFT central charge is positive while the bulk theory is unitary. Of course this works 
only around (anti) de Sitter space and it is not clear how severe the problem with unitarity 
will be around other spacetimes. In any event, it is fair to say that this construction is an 
interesting one and worth probing further. The way that one gets two derivative equations of 
motion in this case is that the six derivative R 3 terms cancel off the offending higher derivative 
terms arising from the four derivative terms. This of course only works around a non-trivial 
background and a similar construction cannot be used around flat space. 

We want to construct an interesting class of higher derivative models in D = 4 such that 
fluctuations around (anti) de Sitter space have two derivative equations of motion. There are 
several motivations behind doing this. 



1. We wish to propose interesting higher derivative lagrangians in D = 4. Had we worked 
just with four derivative R? lagrangians we would be led to the Gauss-Bonnet term 
which does not alter the equations of motion in D = 4 as it is a total derivative. We 
wish to do something more interesting than this. In particular, we want to allow for 
the possibility of exact black hole solutions as in J?]]. However, the approach used in J?J 
does not extend to D = 4. 
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Higher derivative gravity is an interesting playground to consider transport properties 
in interesting field theories at strong coupling . It has been used to study bounds on 
the ratio of shear viscosity to entropy density in the context of quark gluon plasma 



[31, [H|. An unsolved question here is if the ratio can be driven to zero without any 
pathologies on the gravity side. No useful study has been yet carried out in the context 
of d = 2 + 1 CFTs which are useful in the AdS/CMT applications. Part of the reason 
is that no straightforward generalization of Gauss-Bonnet or quasitopological gravity 
exists in D = 4. 

3. We wish to consider c-theorems in this context both in AdS/CFT and cosmology. These 



were studied in higher derivative gravity in [23, 32 1 in the context of AdS/CFT 



4. f(R) theories have been extensively studied as viable alternatives to inflation [[3^, 34]. 



Generalizations to f(GB) or as functions of Gauss-Bonnet have also been considered 



[35|. The f(GB) models may have some problems since they appear to be incompatible 
with observations j33|. Thus (and otherwise) it is interesting to look for alternatives. 

Keeping these motivations in mind we will engineer higher derivative lagrangians includ- 
ing upto six derivative curvature invariants in D = 4. These can be thought to be distant 
cousins of the quasitopological theory of Myers and Robinson |?j . We will find a five parame- 
ter family of such lagrangians which yield two derivative equations of motion for fluctuations 
around (anti) de Sitter spaces and which allow a simple c-theorem as in [^, |2^[. A two pa- 
rameter subspace is found such that the equations of motion for fluctuations around FRW or 
a static domain wall are two derivative. This lagrangian coincides with the choice of Gauss- 
Bonnet for the R 2 terms and a cubic invariant of the Weyl tensor for the i? 3 terms. If one 
considered black holes in such spacetimes, the solutions would receive corrections from the 
six derivative terms. A three parameter subspace is found where exact black hole solutions 
exist. Rather interestingly, we will find that there is a lower bound for the entropy in de 
Sitter spaces. Some aspects of this theory will also be discussed in the context of AdS/CFT 



c-theorems in 32]. 

This paper is organized as follows. In section 2, we write down the equations of motion 
to be used in the rest of the paper and specify our conventions. In section 3, we review 
the construction of new massive gravity and its extensions. We consider c-theorems in the 
context of cosmology. In section 4, we turn to D = 4. After constructing the lagrangian, 
some exact black hole solutions are presented. In section 5, we turn to discussing c-theorems 
in the context of cosmology. We conclude with a discussion of open problems in section 6. 

2. Six derivative theories and equations of motion 

We will be interested in lagrangians that include upto six derivative curvature invariants. 
These take the form M 



2l p 



1 = I d°x^~g 



± {D - l f- 2) + R + L^ + L^ 



(2.1) 
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where 

X4 = XiR abcd R abcd + \2R a bR ab + A3-R 2 , (2-2) 

Z± = fJ>iR a bc d R b edfR aeC ^ + ^RabcdR" 1 ^ R + ^RabcdR^ eR de (2-3) 
HARabcdR aC R M + ^Ra h Rb C Rc a + ^&RaR b R + • 

We have left out terms that involve Vi?'s. There are several reasons for this. Firstly, as shown 
in Q, in order to get exact black hole solutions, these terms should be absent. Secondly, had 
these terms been present, schematically they would look like (VR)(VR) so that when we vary 
this, there would be contributions that look like RW5R so that the equations of motion for 
fluctuations would necessarily involve more than two derivatives unless these contributions 
canceled among themselves. We will assume in what follows that there are no such terms 
although including them should not be a problem so long as their contributions to fluctuations 
cancel among themselves. 

The equations of motion that follow from this are given by J?], [J 

Rab ~ \ 9 abR T {D ~fj^~ 2) 9ab - L 2 K a f - = , (2.4) 

where 

K$ = X 3 (-2RR ab + 2V a V b R + g ab [^R 2 - 2V 2 R\) 

+X 2 (-2R c a R cb + 2V c V (a R c b) - V 2 R ab + g ab [^RcdR cd - \^ 2 R}) 

+^i(-^9abRcdefR cde ^ — 1RacdeR b de ~ 4V 2 i? a b + 2V a V b R + 4R c a R bc + AR cd R c ^ d ) , 

(2.5) 



while 




= M"3^ ec f R C9e aR fg b + SV d V c R e d f{b R a) ecf - 3V c V d R {a e b) f R e d f + ± gab R c d J R d 9 f h R g \ e ) 
+/^2(— 2RR acde R b cde — R ab R cde fR cde f + V '&V a R cde fR cde f + 4V c RR^ d b 
+ l -g a b[RRcdefR cdeS - 2V 2 R cdef R cde f]) 

+fi 3 (-2R acd s R b cdf R e f - R^fcRfoffaRfyc - ^ 7^ ' (bR a ) ecd R c dJ ~ -j^ 2 R cde a Rcdeb 

+ 2V d V c R e %R b)e + 2V d V c R e (ab fR e d + \g a b[R cdef R cde 9 Rf 9 - V f V c R d ^R d J}) 
+H±{—ZR c ( a R b y e R e + 2VdV( 6 i? a ^ c e R ce — V 2 R acbd R c + V d V c R^R b ) — VrfV c i? c R ab 
+ \gab[R cd Rced } R ei - V d V c R e c /R e f]) 

d 3 2 d 1 3 d 

+im(-3RacR C dRb - 2^ RacRb + 3V c V( a i? b )-R^ + g ab [-R c R e d R c e - -V d V c R c e R e ]) 
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+ f i 6 {-R d c R d R ab - 2RR ac R c b + V b V a R cd R cd - V 2 RR ab + 2V c V (b R c a) R 

+ 9ab[\RRcdR cd - V d VcRR cd - V 2 R cd R cd ]) 
+fi 7 (-3R 2 R ab + 3V a V b R 2 + g ab [^R 3 - 3V 2 R 2 ]) . (2.6) 

For computational purposes, it is sometimes easier to use the effective action approach as 
described in jjffl. Suppose one is interested in finding planar black hole solutions. Here one 
starts with an ansatz for the metric of the form 

ds 2 = - N (rff(r)dt 2 + + j^idx 2 + dy 2 ) , (2.7) 

plugs this into the action and works out the equations of motion for N(r) and f(r). When we 
deal with fluctuations around a given background, we turn on all metric components for the 
fluctuations with arbitrary spatial and temporal dependence, expand the action upto second 
order and work out the equations of motion for the fluctuations. 

Finally we make a note of the following relations with the invariants of the Weyl tensor 
Cabcd m D = 4. At quadratic order, there is a unique invariant which is related in the 
following way: 

r-i riabcd tj rtabcd q r> r>ab , ^ r>2 fo Q\ 

^abcd^ — ttabcd-K ~ z ^ab^ + g-K > (,*-°) 

At cubic order there are two invariants 

tj/ /-~irstu/-~ivw s-i tj/ r^rstur^v w s~i in n\ 

VV\ — U U rt^svuw , VV2 — ^ ^ r t^swuv ■ \ z - y ) 

aW\ + (3W2 can be expanded in terms of Z4 by choosing the follows: 

a+J3 a -5/3 11a - 43/3 

Mi = > M2 = — 77— > M4 = -3/3 = -2^3 , /i 5 = 8/i 2 = -Me , M7 = 771 • 

2 16 144 

(2.10) 

Now there is a Schouten identity that leads to the fact that W\ = W2 in D < 5. As a result, 
in order to be consistent, it had better be true that 

v e> rtabcd r> i p r>abc nde i c p pac r>bd ior>br>cr>a o r> b r> a r> , p3 rv 

(2.11) 

Where does this identity come from? This comes from the fact that in 5 dimensions we can 
construct 

e abcde e fghij R ab jr g R cdhi R e j = —4X 5 , (2-12) 

and this should vanish in four dimensions. Thus using this identity we can drop one of the 
/Vs for i > 2. For example, we could choose to drop n^R 3 . In what follows we will set /3 = 
when D = 4 but will retain all the /Vs. This will serve as a cross-check on the algebra. 
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3. Extensions of new massive gravity 



In this section we review the construction of six derivative lagrangians in ||] in D = 3. In 
D = 3, the Riemann tensor is given in terms of the Ricci scalar and Ricci tensor. As a 
result the number of independent cubic terms reduces to just three. These lagrangians were 
constructed so that a simple c-theorem existed on using the null-energy condition. The action 
is given by 



1 

21p~ 



d' i X^ 



-g(R+^ + L 2 TZ 2 + L 4 TZ 3 )^^- 



d 3 x^{R+-^ + K), (3.1) 



where 



Tl 2 = 4(XiR ab R ab + A 2 i? 2 ) , 

TZ 3 = ^(^R b a RlR a c + fi 2 R ab R ab R + ^ 3 R 3 ) . 

Anticipating a relation with the AdS/CFT correspondence, we consider 

ds 2 = e 2A{ - r \-dt 2 + dx 2 ) + dr 2 . 
By demanding that there exists a simple function such that 



c'(r) 



t\ - t; 

' £ P A' 2 



>o, 



it was shown that A2 = — 3/8Ai and /ii = §7/^3,^2 = — ?=/X3. This led to 



c(r) = 

satisfying 

In the absence of a matter sector 



17^d) H'Z — 17/ 



t P A 



7 (l + 2A 1 L^ + / u 3 ^ /4 ), 
c'(r) > 0. 



Air) = ^ 
L 



rfc 



1/2 



L 



1 - foo + /^Ai + /^M3 = 



(3.2) 
(3.3) 

(3.4) 
(3.5) 

(3.6) 
(3.7) 

(3.8) 



That AdS is a solution to a higher-derivative theory is not surprising [37]. After all, the 
only effect (if any) of curvature corrections would be to correct the AdS radius. The four 
derivative theory (fi 3 = 0) has been studied in [19, 20, 21, pq| . Quite curiously, the relative 



coefficients of the R 3 terms work out to be the same as that in the Born-Infeld extension 
considered in [39j. It was further shown in the first paper in [22] that even the R 4 terms work 
out to be the same in the two approaches. Black hole solutions in the six derivative theory 



were considered in |4C]. An infinite order generalization of this construction was shown in 
[pC[| . It was argued in J3(| ||] that by suitably tuning the parameters (in the above example 
choosing Ai = —^sfoo), the equations of motion for fluctuations work out to be second order. 
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In the six-derivative extension considered above, this implies that /oo = 1 or in other words 
the radius of AdS is unaffected. In fact in this case the fluctuations arise from a Fierz-Pauli 



action of the type [42] 

+ f(^V-^ 2 ) + 0(^ 3 )) • (3-9) 

Here V p is with respect to the AdS background metric g^ u with no fluctuations and h^ u 
is a small fluctuation. Here h = h^. In order to see the equivalence, one needs to choose 
the transverse traceless gauge. This is precisely the action one gets when one expands the 
Einstein-Hilbert action around AdS (or de Sitter) in the presence of a cosmological constant 
A (which in our case for AdS is — 1/L 2 ). In other words, the combination of the four and six 
derivative terms are such that the only effect on fluctuations (upto 0(h 2 )) is to renormalize the 
coupling constant. Since we are in D = 3 there are no propagating modes and as such there 
appear to be no restrictions on (1 — ^3). However the CFT central charge is proportional to 
1 — fj,^ as a result we need ^3 < 1 for there to be a sensible AdS/CFT dictionary. Of course, 
it would have been rather weird if this condition was not satisfied since then the effective 
Newton constant would be negative! One other interesting point to note about the above 
construction is that it is crucial to have the cosmological constant to begin with for this to 
work. An easy way to see this is the following: to get rid of the cosmological constant, we need 
to send L — > 00 while rescaling H3L 4 = /x, A1L 2 = A keeping /}, A fixed. Using eq. (|3.8|) this 
would then simply lead to /oo = 0. This seems to be a peculiar feature of the six derivative 



extension and is not expected to hold in general [30|. 

Another point to note in this construction is that if we assume that there is a flow between 
two asymptotically AdS spaces, then the equations of motion for fluctuations are going to be 
two derivative only around one of the spaces. This is because we have tuned the parameters 
to yield two derivative equations of motion for a specific /oo- Once we change /oo, which is 
what will happen in a flow, the equations of motion around the other AdS will no longer be 
two derivative. This is markedly different from what is possible in D = 4 and higher where 
one can tune the parameters to get two derivative fluctuations around any (anti) de Sitter 
space. 

Let us summarize the way the theory was constructed. The motivation above was to have 
a simple c-theorem. Furthermore, we demand that the equations of motion for fluctuations 
around AdS are two derivative. This is the point of view that we will take to construct 
lagrangians in D = 4. If we instead demanded that fluctuations around de Sitter should be 
two derivative, this is also easily achieved. Either one starts anew by flipping the sign of 
the cosmological constant and then studying fluctuations around de Sitter, or one replaces 
L — > iL, r — > it, t — > iy in the AdS solutions. This leads to the constraint Ai = /i3 on the 

parameters as opposed to Ai = —^3 in the AdS case. The form for the action for fluctuations 

2 



is still the same as eq. (|3.9| ) with A = 1/L in this case. 
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One further nice property of NMG (and the above extension) is the following. It was 
pointed out in || that the c-function has a connection with the Wald formula since 

1 BC 

odrrgttj^^- = 1 + 2L 2 X 1 A'(r) 2 + LVs^V) 4 = i P c(r)A'(r) , (3.10) 

1 OK r trt 

where C is the lagrangian. In the de Sitter case, there is a similar relation as well. If we write 
the metric as 

ds 2 = -dt 2 + e^(dx 2 + dy 2 ), (3.11) 

then 

Ig^gtt^— = 1 - 2L 2 \ 1 H(t) 2 + L^ 3 H(t) 4 = ^s(t)H(t) , (3.12) 

2 oRxtxt 2ir z 

where ^{t) = H(t), the Hubble parameter. Using the equations of motion and assuming the 
null energy condition on the additional matter sector, we find that 

s(i)>0. (3.13) 

In other words, s(t) is increasing with time. Now in this case s(t) has the natural interpretation 
of the entropy on the apparent horizon. To see this recall that Wald's formula for entropy 



f dC 

S = -2vr j dx\/h-^—e ab e cd , (3.14) 

where e a & is the binormal to the horizon. In this case S works out to be 

ttA dC 

z> = -^—gxxgttTrs — • (3-15) 

lp Otixtxt 

Now the area of the apparent horizon is simply A = 2ir/H(t) since the location of the apparent 
horizon is at a distance D = \/H{t) and we are in two spatial dimensions. In de Sitter space 
this gives the correct definition of entropy for Einstein gravity [44, 45] on the cosmological 
horizon. This leads to the interpretation of s(t) as the entropy and we see that if the null- 
energy condition is satisfied on the additional matter sector (with respect to the full metric) 
then entropy increases with time. We would want s(t) > 0. First consider /X3 = 0. This 
immediately leads to 

,2.1 ... - n 



or in other words, there is an upper bound for H if Ai > 0. Turning on fi 3 = Ai this changes 
to 



HL<^ Xl + ^V *^l, if A! > 1 . (3.17) 

One curious feature of this theory is that it is not necessary for area to be increasing for 
entropy to be increasing. To see this consider for simplicity [13 = 

s(t) = -j^2H(l + 2L 2 \ 1 H 2 )>0. (3.18) 
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Now in standard cosmology with Ai = /i3 = this implies H{t) < or in other words the 
standard result that the Hubble parameter decreases with time (or A oc 1/H(t) increases 
with time). However in the presence of X\, fJ-3 we see that this no longer holds. If fi3 = then 
with (1 + 2L 2 X 1 H(t) 2 ) < we indeed have H > or in other words the area decreases with 
time but due to the higher derivative contribution which overwhelms this decrease, entropy 
still continues to increase. It will be very interesting to consider this example in the context 
of cosmological bounce solutions 2 . It appears at first sight that this regime will contain ghost 
excitations. 

If we were to consider fluctuations around the standard FRW metric 

ds 2 = -dt 2 + a{tf{dx 2 + dy 2 ) , (3.19) 

then there would be more than two-derivative terms in the equations of motion. All these 
terms come multiplied by either 

(X 1 -L 2 f i 3 H 2 -2L 2 ^H)\, 

gt 



or 

(Ai-LVa^ 2 )^ 



Here H{t) satisfies 

H 2 - A + L 2 X 1 H i - L 4 fi 3 H 6 = 8tt P , H(l + 2L 2 \ l H 2 - 3LV3# 4 ) = -8tt(P + p ) • ( 3 - 20 ) 

Thus in the absence of matter (p = P = 0) if H = y/X = 1/L with Ai = ^3, then the 
potentially problematic three and higher derivative terms would disappear from the equations 
of motion and we would be left with the promised two derivative equations for the fluctuations. 
Furthermore, when we add matter we would expect a(t) ~ t p so that H = p/t in which case the 
ghost terms would be proportional to either (Ai — L 2 ^^ 2 /t 2 )p/t 2p+1 or (Ai + L 2 ^p 2 /t 2 )/t 2p . 
In both cases as long as p > these terms would become irrelevant at large times. Also note 
that as in the AdS case, the fluctuations are two derivative only around one of the de Sitter 
vacua if there is a flow between two de Sitter vacua. 

Let me end this section by pointing out that Schwarzschild de Sitter solutions are straight- 
forward to find. It is easy to check that 

ds 2 = -f(r)dt 2 + i r -+ r 2 dtf , (3.21) 
IV) 

with /(r) = l — 2m—r 2 HQ satisfy the equations of motion where Hq — A+ L 2 XiHq — L 4 = 
0. In de Sitter in general dimensions, there is a maximum mass that the black hole can have 
which is fixed by the fact that in this case the de Sitter horizon and black hole horizon coincide. 



2 The first of eq.(3.2C) suggests that for p > 0, in order to have a bounce A < 0! It is still possible to have a 
de Sitter solution in this case, however the fluctuations will now involve more than two derivatives. Bouncing 
cosmologies in non-relativistic higher derivative theories have been proposed in |36||. 
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In D = 3 there is only the de Sitter horizon When /(r^) = for = \f\ — 2m/Ho, 
there is a conical singularity with a positive deficit angle at the origin corresponding to a 
point-like mass at the south pole. Further there is a maximum value of m = 1/2. The Wald 
entropy is given by 

Swaid = ^-(1 - 2L 2 \ l Hl + LVs^o) , (3-22) 

with A = = 2tt\/1 — 2m /H which further justifies the identification of s(t) in eq.( |3,12 ) 
as the entropy of de Sitter space. 



Comment on the sign of Ai 

Let me briefly comment on the sign of Ai when p^ = 0. There are two interesting possibilities 
for either sign of Ai. If Ai > 0, then (1 + 2L 2 \ 1 H 2 ) > and we have the sign of the kinetic 
term to be positive and an upper bound on the Hubble parameter set by eq.( p,16 ) arising 



from demanding positive entropy. If Ai < then (1 + 2L 2 XiH 2 ) will run the risk of turning 
negative as H increases and introducing ghosts. In this case there will be no upper bound 
on the Hubble parameter although entropy is positive. Furthermore combined with eq.( |3.16 ) 



or eq. (|3l7D and the first of eq. flpg ) it is easy to see that a(t) cannot be zero for radiation 



or dust matter. In this case p ~ 1/a 2 or p ~ 1/a 3 which blows up. Thus to be consistent 
H — > oo but this is forbidden by an upper bound on H{t). As such aesthetically it may seem 
that Ai > is a more pleasing possibility! It will be interesting to understand what goes 
wrong when the upper bound on H(t) does not hold and s(t) becomes negative. 

4. Higher derivative gravity in D = 4 

It is interesting to ask the following: Given an action for the fluctuations around (anti) de- 



Sitter of the form eq^^), what is the covariant action whose quadratic expansion, 0{h 2 ) 
would give rise to this? If we worked using a truncated lagrangian which included upto only 
four derivative terms, the answer is that only Einstein-Hilbert with a cosmological constant 
would work. This is because R 2 would generally lead to four derivative equations of motion. 
In D = 3 this is not necessarily problematic: the new massive gravity is consistent under 
certain circumstances. However, in higher dimensions only Gauss-Bonnet or Lovelock terms 
would lead to ghost-free theories. In D = 4 this is not interesting since Gauss-Bonnet would 
not contribute to the equations of motion (however, the Gauss-Bonnet terms can contribute 
to the entropy in a crucial way as we will discuss in the next section). Although f(GB) 



theories have been considered in the literature (see [33] for references), these theories are not 
viable for cosmology. This gives us motivation for searching for other creative theories. 

In D = 4, if one expanded around flat space, the equations of motion for a generic four 
derivative gravity theory will always be more than two derivatives 3 . It is in fact easy to argue 



3 For R + R 2 , the theory can be mapped to Einstein gravity and scalar field which is free of the usual 
problems |B3|. 
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that adding even higher curvature terms would be of little use since the equations of motion 
would typically involve hR 2 ~ 0. However, around a background like (anti) de Sitter, the 
story is not so simple. In fact as we saw in the D = 3 example, one can play the R 2 terms 
off the i? 3 terms to have two derivative equations of motion. We are now going to show that 
there is a similar construction in D = 4. One nice feature of this theory, unlike the D = 3 
example, is that one can have two derivative fluctuations around any (anti) de Sitter space. 



Our starting point is the action given in eq.(2.1) in D = 4. Although we have left an 
explicit cosmological constant in the lagrangian, we will soon argue that this is not necessary 
to have a (anti) de Sitter solution with two derivative fluctuations unlike the D = 3 case. Let 
us begin with the AdS case. The metric is 

ds 2 = e 2A( * r \-dt 2 + dx 2 + dy 2 ) + dr 2 , (4.1) 

where A(r) = r\J foo/L. The parameters need to satisfy the following constraints in order 
for the equations of motion for the fluctuations to be two derivative 4 and such that a simple 
c-theorem exists: 

Ai + A 2 + 3A 3 = 0, (4.2) 
3/ii + 48/x 2 + Ufj, 3 + 16^4 + 18^5 + 60^ 6 + 216/i 7 = , 
12/i 2 + 4/i 3 + 3/x 4 + 5/i 5 + 12/i 6 + 36^7 = , 
A 2 + 4Ai - 3/oo(-8^ 2 - 2/i 3 - fj,4 + /U 5 ) = . 

Note that using eq.( |2.11| ) we could set fij = 0. The first three constraints will be imposed 
explicitly while we will leave the last one implicit in what follows. By imposing all the 
constraints we will find theories that have second order linearized equations of motion and 
admit a simple c-theorem. The first three constraints were obtained by demanding that a 
simple c-theorem exists as in the D = 3 case. The last constraint is needed to make the 
linearized equations of motion two derivative. This will also serve the purpose of illustrating 
that it is not necessary to demand two derivative equations of motion for fluctuations to 
have c-theorems. The first of the constraints is satisfied both by Gauss-Bonnet and by Weyl- 
squared. In the absence of the six-derivative terms, the last constraint would have uniquely 
selected Gauss-Bonnet as one would have naively guessed. However, the interesting thing to 
note here is that in the presence of six-derivative terms, it is not necessary for the R 2 terms 
to be Gauss-Bonnet to have two derivative fluctuations. Of course, this is a highly contrived 
situation where the couplings Aj,/Vs have to be fine-tuned to make things work. If the R 2 
and R 3 couplings were independent of one another then we would be forced to choose the 
Gauss-Bonnet term. Furthermore, using Gauss-Bonnet and W± as defined in eq.( |2.9| ) will 
satisfy all the constraints as well. If we were to use W\ however, = — AL 2 /3 and there 



4 In a previous version of this paper it was erroneously implied that these constraints are necessary for 
two derivative equations of motion for the fluctuations. While these are sufficient conditions, these are not 
necessary. In fact it is shown in 0], that it is possible to have two derivative equations of motion without 
having a simple c-theorem. 
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is no way to get AdS without an explicit cosmological constant in the first place. We would 
like to keep things as general as possible, not only because of this, but also to explore the 
possibility of having exact black hole solutions. Another point to note is that with W\ the 
c- function would not receive corrections from the six derivative terms. So we will keep things 
more general for now. Here in the absence of matter /oo satisfies 

l-/oc-(w + M4-M5)/^ = 0. (4.3) 

Thus we have a five-parameter family of theories. Using the equations of motion and the null 
energy condition one can show that c(r) defined through |}2| 

1 dC 

^grrgttjTTT- = [1 + 2(2Ai + \ 2 )L 2 A' 2 - 3( W + ^ - H5)L A A' 4 )] = £ 2 P A' 2 c(r) , (4.4) 

is monotonically increasing. In other words cjjy > cm when there are UV and IR fixed 
points. Another point to note is that as in 0, one can infer unitarity from the left hand 
side of eq.( |4.3| ). The action for the fluctuations take the form (1 + 3(/ii + /i4 — ^)foo 2 ) times 
what we have in Einstein gravity. As a result, (1 + 3(//i + \i± — ^)foo 2 ) > needs to hold for 
unitarity. Together with eq.( [4.3D this translates into the condition fi% + ^4 — ^5 > —4/27. 

As in the D = 3 case one will have two derivative fluctuations only around one of the 
AdS vacua if there is a flow between two different AdS vacua. However unlike the D = 3 case, 
we have a bit more freedom here. In fact in addition to eq.(|4.2|) if ^5 — ^4 = 8fi 2 + 2yU3, then 
the fluctuations around any AdS vacua will be two derivative! In this case, the last constraint 
in eq. fl4.2|) gives us the familiar result that the R 2 combination will be Gauss-Bonnet (and 
hence topological). 

In fact we can do much better. It is possible to actually have the fluctuations around any 



A(r) in eq.(4.1) to be two derivative and in fact the same as in Einstein gravity. In a sense 
this is exactly like Gauss-Bonnet but at six derivative order! The choice of parameters that 
accomplishes this is 

A 2 = -4Ai = -4A 3 , fi 7 = ^ + H2, M5 = 8/J.2 = -pe, ^4 = — A*i + 8/i 2 = -2^ 3 • (4.5) 

3d 



If we use eq.(2.11) and set ^7 = then we will have a two parameter family of lagrangians 
with nice properties. Now, unlike Gauss-Bonnet which is topological and does not enter any 
equation of motion in D = 4, the six derivative terms are more interesting. Only under 
further conditions does one find simple exact black hole solutions as we will mention shortly. 
In other words, the six derivative terms are not topological. With this choice of parameters 
in empty AdS, /oo = 1. This finding will continue to hold for fluctuations around de Sitter 
as well. Also note that aW\ satisfies the constraints in eq.( [4.5[) . However, in what follows, 
unless explicitly specified, we will only impose the first three constraints in eq. fl4.2|) . 

In order to see that it is not necessary to have the explicit cosmological constant in the 
theory, we rescale Aj = AjL 2 ,/ij = ^L 4 and send L — > 00. Now the AdS radius works out to 
be foo/L 2 = /oo. In this case, foo satisfies 

/oo + (Al + A4-A5)/i = 0. (4.6) 
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Now it is possible to have foo = 1/ \/7i5 — /"4 — Ai as a non-trivial solution. However, as in 
Boulware and Deser j5| fluctuations around these vacua will contain ghosts. Note that 
does not depend on the R 2 parameters although the fluctuations do. 

It is also possible to get exact black hole solutions 5 if in addition to the first three 
constraints in eq.( |4.2D we restrict /Ui = 0, 4/i 2 + /i3 = 0. In this case it is possible to show that 

^7^2 ~.2 

ds 2 = -f{r)dt 2 + ^- + — 2 {dx 2 + dy 2 ) , (4.7) 
J[r) L z 



with f(r) = fooT 2 /L 2 (l — M/r 3 ), with obeying eq.( [4.3D , satisfies the equations of mo- 
tion. The equations of motion involving f(r) are four derivative. In spite of this we have 
a remarkably simple solution! In fact apart from the corrected asymptotic AdS radius, the 
solution is of the same form as in Einstein gravity. We do not yet have an understanding for 
this simplicity. The complicated nature of the equations of motion also makes the proof of 
Birkhoff 's theorem along the lines of |J a difficult problem. We have a three parameter family 
of exact solutions. It is not possible to use Weyl invariants to accomplish this since setting 
Hi = would need a = —f3 in aW\ + fiWi which vanishes as it is a Schouten identity. We 
have here more general exact black hole solutions where /oo 7^ 1 (if we use the last constraint 
in eq.(|4,2|) then = 1). The ratio of shear viscosity to entropy density rj/s in this model 
can be calculated following for example [jlCR . It works out to be 1/4-7T in spite of the fact 
that there are higher derivative terms, rj and s each get corrected but in the same manner. 
It will be interesting to constrain the parameter space following and figure out which 
hydrodynamic quantities do get corrected. We will leave the possibility of exact black holes 
in theories with /xi ^ 0,//3 7^ — 4// 2 as an interesting open question. 
We can easily extend the analysis to de Sitter space. We find that 

ds 2 = e 2H °\dx 2 + dy 2 + dz 2 ) - dt 2 , (4.8) 

is a solution with Hq satisfying 

AL 2 



H$L Z - (/ii + fi 4 - V5)H$L b = . (4.9) 

To get two derivative equations of motion for the fluctuations, the parameters satisfy the 
following constraints: 

Ai + A 2 + 3A 3 = 0, (4.10) 
3/ii + 48/x 2 + 14/u 3 + 16^4 + I8/X5 + 60/u 6 + 216/x 7 = , 
12/i 2 + 4/i 3 + 3// 4 + 5/i 5 + 12/i 6 + 36/i 7 = , 
A 2 + 4Ai + 3F 2 L 2 (-8^ 2 - 2/i 3 - M + Ms) = . 



5 These solutions would not have been found using the method in Q since the equations of motion involving 
f(r) in our case is fourth order but an exact solution nonetheless exists. The method in Q will only work 
when the /(r) has second order equations of motion. 
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We will impose the first three of these constraints explicitly and leave the last one implicit. 
Furthermore, if ^5 — /U4 = 8^2 + 2/i3, then the fluctuations around any de Sitter will be two 
derivative while if eq.(i.E) holds then the fluctuations around any FRW will be two derivative. 
As before we will not impose these constraints unless explicitly specified. Again we can rescale 
the original cosmological parameter out by taking L — > 00. Now Hq = 1 / yjjl^ — jl^ — fi\ 
although the resulting vacua will again contain ghosts. One can also construct exact asymp- 
totically de Sitter Schwarzschild solutions if [i\ = 0, 4/^2 + = 0. In this case it is possible 
to show that 

dr 2 



ds 2 



-f(r)dt z + 



f(r] 



+ r 2 (d9 2 + sin 2 



(4.11) 



is an exact solution with f(r) = (l—2m/r—r 2 H 2 ) with Hq satisfying eq.((4.9|). The asymptotic 
space, r — > 00 is de Sitter in static coordinates |4q| . The Wald entropy 6 of these black holes 
works out to be 



S- 



Wald 



2ttA 


4L 2 m\i 


£% 


r h 


2ttA 


AL 2 m\i 


f 2 


r h 



+ {1 - 2(2Ai + \ 2 )H 2 L 2 - 3(^4 - ^)H^L 4 } 



+ s £ P H 2 



(4.12) 
(4.13) 



where A = 4-7rr 2 is the area of the horizon and is the location of the black hole horizon. 
The second equation defines sq whose meaning will be clear in the next section. Now de Sitter 
space has an interesting property. There is a maximum mass black hole that one can fit inside 
de Sitter [^]. This happens when the de Sitter horizon and the black hole horizon coincide. 
The location of the horizon is found by solving the depressed cubic equation f(r) = 0. When 
m = 1/(3\/3Hq), then there are two equal positive roots = 1/(\/3Hq), in other words the 
de Sitter horizon and the black hole horizon are on top of each other 7 . Thus demanding that 
the black hole entropy is positive we get the following interesting inequality 



4L\ 
so > -p- Ai . 



(4.14) 



In the next section we will see that so is the entropy of de Sitter space. Thus we seem to 
find that there is some minimum entropy in the system if Ai > 0! In D = 3, there are 
exact Schwarzschild de Sitter solutions as well (with no constraints on Ai or ^3) but their 
Wald entropy does not receive any contribution analogous to the \\m term in D = 4. Hence 
eq.(4.14) seems special to D = 4. 



5. Entropy theorems in cosmology 

Since c-theorems in the context of AdS/CFT will be discussed extensively in (§3], we will focus 
on c-theorems in the context of cosmology. But before we begin, let me briefly summarize 

6 The fact that Gauss-Bonnet contributes to Wald entropy in 4 dimensions was pointed out in a different 
context in 

7 m/rf l is maximized for this choice as well. 
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the findings in [p2f| . In |23|], we found that in the context of the quasitopological gravity 
introduced in 0, there are nice c-theorems in arbitrary dimensions. In odd dimensional 
CFTs (D = 4 for e.g.) the quantity that is flowing was interpreted as entanglement entropy 
of the CFT on S D ~ 2 x R. In [32], it is shown that these c-theorems can be extended for a 



wide class of higher derivative gravity theories. In the context of cosmology, we naturally 
interpret the quantity which flows as the entanglement entropy in de Sitter space across the 
cosmological horizon. Let us consider this calculation using the static coordinates since in 



these coordinates, the results can be compared directly with those presented in [23] in the 



AdS/CFT context. The calculation can be done using the metric eq.( 4.1l| ) which leads to 



eq. (|4.12[ ). The de Sitter entropy is extracted by setting m = in eq.( [1.12 ). The interpretation 



of the AdS/CFT calculation was that the c- function at the fixed points gives information 
about the entanglement entropy between the two halves of S 2 when the CFT is placed on 
S 2 x R. In the de Sitter context, the interpretation 8 on the bulk side is straightforward. It is 
simply the entanglement entropy of the two disconnected regions 9 of de Sitter space. 

Now let us turn to c-theorems in cosmology. For this we want to start with a metric of 
the form 

ds 2 = -dt 2 + e^ t] {dx 2 + dy 2 + dz 2 ) . (5.1) 

Now we calculate 

1 dL I 2 

^gxxgttj^ = 1 - 2(2Ai + X 2 )L 2 H(t) 2 - 30: + ^ 4 - ^)L A H{tf = -JL s (t)H(t) 2 , (5.2) 

2 oRxtxt 4vH 

where s(t) is a monotonic function which follows from a straightforward application of the 
equations of motion and assuming the null energy condition (assuming H{t) > 0). We see 
that this is in agreement with eq.( 4.12j ) in the zero mass limit as expected. One important 



thing to note here is the following. Let us impose eq. (|4.5| ) so that the six derivative terms do 
not contribute to the entropy and the R 2 combination is just the Gauss-Bonnet term. The 
equations of motion do not get affected due to this term. However, the entropy does. The 
entropy here is given by 

4tt 2 

s(t)\ G B = ^p^p(! + 4A 1 L 2 tf(i) 2 ) . (5.3) 

As in the new massive gravity example, we have an upper bound for H < ^J-l/^) if 
Ai < 0. It will be interesting to ask what goes wrong when the entropy becomes negative. 
The analysis in [46] may have some relevance in this case. Note that eq. Q5.3j ) also implies 

s(t)\ G B ■ 

4tt 2 



sq > 4L 2 Ai/^p which is exactly what we found in eq.( [4.14 ). Entropy of de Sitter 



spaces have been investigated previously in four derivative theories in general dimensions in 



following [47] and our conclusion is in agreement with their analysis. 



8 In the dS/CFT context, it is again mapped onto the entanglement entropy of the two halves of S 2 . 
9 The horizons are at r = 1/Hq. As discussed in 0], no single observer can access the entire de Sitter 
spacetime and the entropy in this context is the entropy of whatever is behind the horizon. 
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6. Discussion 



We considered higher derivative lagrangians in D = 3 and D = 4 which were constructed such 
that fluctuations around (anti) de Sitter had two derivative equations of motion and such that 
these theories admit simple c-theorems both in the context of AdS / CFT and cosmology along 
the lines of || |23|. |3^]. The most general such lagrangian in D = 4 [see eq,(|2.1|) subject to 
eq. Q4.2D or eq.( 4.10| )] had five parameters. In addition if we imposed ^5 — ^4 = 8//2 + 2//3, then 



the fluctuations around any (anti) de Sitter space have two derivative equations of motion. 
Exact black hole solutions 10 were found if \i\ = 0, 4/i 2 +^3 = along with eq.Q4.2p or eq. fl4.10p . 
Unlike Q these black holes were similar what to what is found in Einstein gravity except that 
the asymptotic (anti) de Sitter radius gets corrected. It will be interesting to analyze these 



theories in detail as in [10 1. As was noted in (7[] ; the construction outlined there does not 



extend to D = 4 or D = 6. The approach used in this paper following [32| is a useful way to 
extend their theories to these dimensions. In the course of our work, we also found that if we 
asked for two derivative equations of motion for fluctuations around any FRW background or 
around a static domain wall in D = 4, then the most general lagrangian upto six derivative 
curvature invariants is L = y/—g{R — 2 A + XGB + aW\) where W\ is defined in eq.( |2.9[ ). 
There is no analogous solution 11 in D = 3. 

In D = 4 an inequality was found which suggests that there is a lower bound for the de 
Sitter entropy. This arises due to the contribution of the R a }, c dR abcd term in the lagrangian and 
is present even though the equations of motion for the background do not receive contributions 
from this term. This is a curious feature of D = 4 and it will be interesting to investigate an 
analogous feature in D = 6 due to R? terms which also will not contribute to the equations of 
motion but will contribute to the entropy via Wald's formula. Of course fluctuations around 
the background may be sensitive to these terms in general except when the combination is 
Gauss-Bonnet in D = 4 or Lovelock in D = 6. The lower bound is sensitive to the sign of 
the R a i, c dR abcd term and exists only if the coefficient is positive. In string theory, scattering 
amplitude calculations in heterotic or higher curvature corrections to the D-brane effective 
action both lead to positive coefficients f|, 52]. This makes a closer analysis of this bound 
very interesting. At the level of this paper, it is not clear how general this bound is or if it is 
specific to the set of lagrangians permitting exact black holes considered in here 12 . 

It will also be interesting to carry out a more systematic study of perturbations using 
these lagrangians. If eq.(EO) is not imposed, perturbations around FRW have fourth order 



10 Actually exact black hole solutions of the type discussed here exist even if eq.(4.2) or eq.(4.1C) were not 
imposed and only pi = 0, 4^2 + M3 = were imposed. 

"In D — 3 the role of the Weyl tensor is played by the Cotton tensor C ab = € acd V c (R b d - \& b d R). Had 
we considered theories made of invariants of the Cotton tensor alone, then it is possible that the fluctuations 
around a general FRW would be two derivative. However, there is no way of making the R 2 contributions to 
the fluctuations cancel on their own. 

12 In the context of AdS/CFT, the effect of topological terms has been investigated in [^, In 4 dimen- 
sions, it was found that the Gauss-Bonnet coupling is fixed to a specific (positive) value. The same finding 
holds for Lovelock terms in higher dimensions. 
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equations of motion, in general we expect there to be ghost modes for these lagrangians. 
However, as we explicitly showed for D = 3 the coefficients of these higher order terms 
become more and more irrelevant as the universe expands. As such, there is a possibility that 
these many of these models may give sensible cosmology. 

We are taking a bottom-up viewpoint |7|, |l0[ for the construction of these higher derivative 
lagrangians as there is very little to go by from fundamental theories like string theory. The 
rules that sensible lagrangians must obey are not clear. For instance, imposing the null 
energy condition is an assumption which need not be true. However, seeing that the entropy 
of de Sitter satisfies a monotonicity property when this holds, makes it worthwhile to think 
that this is a feature that sensible lagrangians must obey. Furthermore, demanding that the 
fluctuations satisfy two derivative equations of motion is certainly not a prerequisite for the 
c-theorems. In fact these theorems hold even though the last equations of (|4.2| ) and ( [4.10 ) 
are not obeyed. As such it seems like a worthwhile pursuit to investigate these models further 
with the hope to learn how to constrain them. 
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